Adhesion between elastic solids with randomly rough surfaces: comparison of 
analytical theory with molecular dynamics simulations 
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The adhesive contact between elastic solids with randomly rough, self affine fractal surfaces is 
studied by molecular dynamics (MD) simulations. The interfacial binding energy obtained from the 
simulations of nominally flat and curved surfaces is compared with the predictions of the contact 
mechanics theory by Persson. Theoretical and simulation results agree rather well, and most of 
the differences observed can be attributed to finite size effects and to the long-range nature of the 
interaction between the atoms in the block and the substrate in the MD model, as compared to 
the analytical theory which is for an infinite system with interfacial contact interaction. For curved 
surfaces (JKR-type of problem) the effective interfacial energy exhibit a weak hysteresis which may 
be due to the influence of local irreversible detachment processes in the vicinity of the opening crack 
tip during pull-off. 



Surface forces play an important role in modern tech- 
nology which deals with micro/nano scale devices. Elec- 
tromechanical devices on this length scale have under- 
gone a rapid development within the last decade. The 
predominant effect of surface forces is due to the increase 
of the ratio between the number of atoms on the surface 
and that in the volume, as the size of an object decreases. 
When two surfaces are brought together, attractive (or 
repulsive) forces act between them, and a non-zero force 
is often required to separate two solid bodies placed in 
intimate contact |l|-l7j , a phenomenon referred to as ad- 
hesion. 

Adhesion manifest itself in different ways. On one hand 
it makes it possible for a Gecko to walk on the ceilings 
or run on a vertical wallJH, On the other hand, adhe- 
sion can lead to the failure of micro or nano devices, e.g. 
micro-sized cantilever beams [10]. Thus, if it is too long 
or too thin, the free energy minimum state corresponds 
to the cantilever beam (partly) bound to the substrate, 
which leads to the failure of the device. However, if the 
surface roughness is increased, the non-bonded cantilever 
state may be stabilized due to the decrease of cantilever- 
substrate binding energy. 

Even the weakest force of interest in condensed mat- 
ter physics, namely the van der Waals force, is relative 
strong on a macroscopic scale. Thus for two bodies with 
perfectly flat surfaces in contact over 1 cm 2 , a force equiv- 
alent to the weight of a car would be necessary to sep- 
arate the surfaces even if only the weak van der Waals 
interaction act between the surfaces. In practice this is 
not the case, and this is referred to as the adhesion para- 
dox. Thus the fundamental problem is not why adhesion 
sometimes is observed, but rather why it usually does not 
manifest itself in everyday life. 

The explanation for the adhesion paradox is that no 
surface of practical use is perfectly flat but have surface 
roughness on many different length scale. In this case the 
interfacial bond breaking will not occur uniformly over 
the contacting interface, but will start at some defect 
and spread by interfacial crack propagation. The high 



stress concentration at the tip of cracks result in a much 
smaller force necessary for the separation of the surfaces 
than would be expected in the ideal case of perfectly 
smooth surfaces and uniform bond breaking. 

In a classical study, Fuller and Tabor |6j have studied 
the influence of surface roughness on adhesion. They 
used silicon rubber balls in contact with (Plexiglas) sur- 
faces with different amount of surface roughness pro- 
duced by sand blasting. They found that with increas- 
ing surface roughness, the pull-off force dropped rapidly 
to zero, and using elastically softer rubber resulted in 
stronger adhesion than for stiffer rubber. This experi- 
mental work, and other similar studies [7], have stimu- 
lated much work to understand the role of roughness on 
adhesion. 

Many practical application of adhesion involves soft 
elastic solids, e.g., pressure sensitive adhesives usually 
consist of weakly crosslinked rubber, which may exhibit 
very complex processes during pull-off such as cavitation 
and stringing [ill [l2|. Here we consider the influence of 
surface roughness on adhesion for the most simple case 
of a purely elastic solid in contact with a hard randomly 
rough surface. If this simplest case cannot be understood 
in detail there is no hope to understand much more com- 
plex systems involving complex rheological materials in 
contact with randomly rough surfaces. 

In order for two elastic solids with rough surfaces to 
make adhesive contact it is necessary to deform the sur- 
faces elastically, otherwise they would only make con- 
tact in three points and the adhesion would vanish, at 
least if the spatial extent of the adhesion forces is ne- 
glected. Deforming the surfaces to increase the contact 
area A result in some interfacial bonding — A7A (where 
A7 = 71+72 — 712 is the change in the interfacial energy 
per unit area upon contact), but it cost elastic deforma- 
tion energy U e \ which will reduce the effective binding. 
That is, during the removal of the block from the sub- 
strate the elastic compression energy stored at the inter- 
face is given back and helps to break the adhesive bonds 
in the area of real contact. Most macroscopic solids does 
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not adhere with any measurable force, which imply that 
the total interfacial energy — + U e \ vanish, or nearly 
vanish in most cases. However, not all the stored elastic 
energy U e \ may be used to break adhesive bonds during 
pull-off but some fraction of it may be radiated as elastic 
waves (phonons) into the solids. This would result in an 
increase in the effective interfacial binding energy, and 
would result in adhesion hysteresis. 

In this letter we study the variation of the effective in- 
terfacial binding energy with the surface roughness am- 
plitude, for elastic solids with randomly rough surfaces. 
We assume perfect (linear) elasticity and compare the 
results of MD-simulations for nominally flat and curved 
surfaces with a recently developed contact mechanics 
theory pjjfi3|. The theoretical and simulation results 
agree rather well, and most of the differences observed 
can be attributed to finite size effects and to the long- 
range nature of the interaction between the atoms in the 
block and the substrate in the MD model, as compared 
to the analytical theory which is for an infinite system 
with interfacial contact interaction. For curved surfaces 
( JKR-type of problem) the effective interfacial energy ex- 
hibit a weak hysteresis which may be due to the influence 
of local irreversible detachment processes. This study 
represents the first test of the theory prediction for the 
effective adhesion energy for 3D systems [17]. 

We review the contact mechanics theory of Persson 
briefly. It can be used to calculate the stress distribution 
at the interface, the area of real contact and the average 
interfacial separation between the solid walls 13|, LL4|. In 
this theory, the interface is studied at different magnifi- 
cations = L/X where L is the linear size of the system 
and A the resolution. The wavevectors are defined as 
q = 2tt/X and qL = 2tt/L so that £ = q/qL- 

Consider an elastic block with a fiat surface in adhesive 
contact with a hard substrate with a randomly rough 
surface. Let <r(x, Q denote the (fluctuating) stress at the 
interface between the solids when the system is studied 
at the magnification £. The distribution of interfacial 
stress 



P(a,0 = (S(a-a( X ,0))- 



(1) 



In this definition we do not include the J(<r)-contribution 
from the non-contact area. 

For perfect (or complete) contact it is easy to show 
that P(cr,C) satisfies 



where 



dP d 2 P 
— = f(C)—— 



f(() = ^E* 2 q L q 3 C(q). 
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Here E* = E/(l — v 2 ) is the effective elastic modulus. 
The surface roughness power spectrum 
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where z — h(x) is the surface height at the point x = 
(x,y) and where (..) stands for ensemble average. The 
basic idea is now to assume that (2) holds locally also for 
incomplete contact. 

To solve (2) one needs boundary conditions. If we 
assume that, when studying the system at the lowest 
magnification £ = 1 (where no surface roughness can 
be observed, i.e., the surfaces appear perfectly smooth), 
the stress at the interface is constant and equal to p = 
Fn/Ao, where is the load and Ao the nominal contact 
area, then P(a, 1) = 5 (a — p). In addition to this "ini- 
tial condition" we need two boundary conditions along 
the cr-axis. Since there can be no infinitely large stress 
at the interface we require P(cr,() — » as a — » 00. For 
adhesive contact, which interests us here, tensile stress 
occurs at the interface close to the boundary lines of the 
contact regions. In this case we have the boundary con- 
dition P{— cr a ,C) = 0, where a a > is the largest tensile 
stress possible. The detachment stress &&(() depends on 
the magnification and can be related to the effective in- 
terfacial energy (per unit area) 7 e ff(C) using the theory 
of cracks fl6| 



leff (C)Eq \ 1/2 



where 



7 eff (C)^*(C) = A 7 A*(Ci)-E/ el (C), 



(3) 



where A*(() denotes the total contact area at the mag- 
nification £, which is larger than the projected contact 
area A(Q. U e \(() is the elastic energy stored at the in- 
terface due to the elastic deformation of the solids on 
length scale shorter than A = L/C? necessary in order to 
bring the solids into adhesive contact (see below). 

From (1) it follows that the area of apparent contact 
(projected on the x?/-plane) at the magnification £, A(Q, 
normalized by the nominal contact area Ao, can be ob- 
tained from 
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da P(a, C). 



(4) 



We denote A(Q/Ao = P P (q), where the index p indicates 
that A(()/Ao depends on the applied squeezing pressure 
p. The area of (apparent) contact at the highest magnifi- 
cation C = Ci gives the real contact area. For the elastic 
energy U e \ we usejlq 

U el « A E^ f * dq q 2 W(q,p)C(q), (5) 

where qL and q\ are the smallest and the largest surface 
roughness wave vectors, and[18| 

W(q,p)=P p (q) [/?+(l-/3)p2( g )] , 
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where /? = 0.45. The equations given above are solved as 
described in Ref. [15j. 

Let us provide some details about the numerical sim- 
ulations. The molecular dynamics system has lateral di- 
mension L x = N x a and L y = N y a^ where a is the lattice 
spacing of the block. In order to accurately study con- 
tact mechanics between elastic solids, it is necessary to 
consider that the thickness of the block is (at least) of 
the same order of the lateral size of the longest wave- 
length roughness on the substrate. We have developed a 
multiscale MD approach to study contact mechanics [l9|. 
Periodic boundary condition has been used in xy plane. 
For the block N x = N y = 512, while the lattice space 
of the substrate b = a. The mass of the block atoms 
is 197 a.m.u. and the a = 2.6 A. The elastic modulus 
and Poisson ratio of the block are E = 77.2 GPa and 
v — 0.42. For self-affme fractal surfaces, the power spec- 
trum has power-law behavior C(q) ~ q~ 2 ^ H+1 \ where 
the Hurst exponent H is related to the fractal dimen- 
sion Df of the surface via H = 3 — Df . For real surfaces 
this relation holds only for a finite wave vector region 
qo < q < q\. Note that in many cases, there is a roll- 
off wave vector qo below which C(q) is approximately 
constant. Here qL = 27r/L,go = 3qL,qi = 512^. The 
physical meaning is that by choosing q = 3qL one can 
obtain a self-average equivalent to an average over 9 in- 
dependent samples. In MD simulations, the substrate 
is rigid and fractal with fractal dimension Df = 2.2 and 
root-mean-square roughness is varied from h rms = 0.1 nm 
to 2.5 nm. We also studied adhesion where the substrate 
was curved into a nominally spherical cup with the radius 
of curvature R = 4410 A. The calculations are carried 
out at temperature T = K. 
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FIG. 1: The surface topography of a randomly rough surface 
with the root-mean-square roughness 0.3 nm and the Hurst 
exponent H — 0.8 (i.e., fractal dimension Df — 2.2). 

The atoms at the interface between block and substrate 
interact with the potential 



V(r) = 4e 



where r is the distance between the pair of atoms. The 
parameter e is the binding energy between two atoms at 
separation r = 2 1 / 6 ro- In the calculations presented be- 
low we have used the ro = 3.28 A and e = 100.0 meV. 
In order to calculate the interfacial binding energy, first 
we move the (upper surface) block at a low constant ve- 
locity towards substrate until the total force on the bot- 
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FIG. 2: (a) The interfacial binding energy (in units of 
the interfacial binding energy A7 = 71+72—712 for flat sur- 
faces) as a function of the root-mean-square roughness. The 
blue and green curves are the result of the MD-simulations 
obtained during contact formation (approach) and contact 
breaking (pull-off). The red curve is the theory prediction (us- 
ing elastic continuum mechanics) for an infinite system with 
interfacial contact interaction, (b) The area of real contact 
A as a function of the root-mean-square roughness calculated 
using the contact mechanics theory. The substrate surfaces 
are self-affine fractal (randomly rough) with the Hurst expo- 
nent H = 0.8. 



torn layer atoms of the block vanishes. At this point the 
work done by the substrate on the bottom layer of the 
block is defined as j e ffAo. In a similar way we calculate 
7 e ff during the separation of the block from its equilib- 
rium separation to infinite separation (pull-off). We find, 
as expected, that to within the accuracy of our calcula- 
tions, 7eff obtained on approach and during separation 
are identical. For a flat substrate A7 turns out to be 
w 100 meV/A 2 . 

In Fig. [2ja) we show the interfacial binding energy 
7 e ff, in units of the interfacial binding energy A7 = 
7i + 72 — 712 for flat surfaces, as a function of the root- 
mean-square roughness. The blue and green curves are 
the result of the MD-simulations obtained during con- 
tact formation (approach) and contact breaking (pull- 
off). The red curve is the theory prediction (using elastic 
continuum mechanics) for an infinite system with wall- 
wall contact interaction. The general shape of the curves 
are the same, but in the MD calculation the adhesion 
extend to larger surface roughness which we attribute to 
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the small system size and the long range wall- wall inter- 
action used in the MD model. That is, in the MD model 
there is an attractive force even between the non-contact 
surfaces which is absent in the continuum mechanics the- 
ory. This long-range attractive interaction is particular 
important for small systems and for surfaces with small 
amplitude roughness. 

We note that the question of interaction energy at a 
distance is related to the difference between JKR [5] and 
DMT theories [20]. However, the situation is more com- 
plex in our case because here we have surface roughness 
which may generate "large" surface areas where the solid 
walls are closely spaced. In the JKR and DMT theories 
(involving smooth curved surfaces) the surfaces are only 
closely spaced close to the rim of the (circular) contact 
area. In our continuum mechanics theory we could in 
principle take (approximately) into account the interac- 
tion between the non-contact surface area using the (cal- 
culated) distribution of interfacial separations P(ix)[2l(. 

In Fig. Efb) we show the calculated area of real contact 
A, at the point where the external load vanish, as a func- 
tion of the root-mean-square roughness calculated using 
the contact mechanics theory. Note that, as expected, A 
vanish for the same roughness amplitude where 7 e ff = 0. 

In an earlier paper [7[ it has been shown that the pull- 
off force between silicon rubber balls and surfaces with 
different types of (non-fractal) roughness could be ex- 
plained by a very simple theory where 7 e ff was calculated 
assuming that complete contact (A = Aq) occur between 
the rubber and the substrate in the nominal (or appar- 
ent) contact area. Only for large roughness, where the 
pull-off force was ~ 20% or less of the pull-off for the 
flat substrate, did this approach underestimate the pull- 
off force. This is consistent with the theory presented 
above. Thus in Fig. [3] we show the interfacial binding 
energy as a function of the root-mean-square rough- 
ness using the full theory (red line, from Fig. [2]), and the 
result obtained assuming full contact for all rms rough- 
ness values (green line). Only for j e ff/Aj < 0.2 does the 
full-contact theory underestimate the interfacial binding 
energy, which is consistent with the experimental data 
presented in Ref. [7[. 

Most solid objects have some macroscopic curvature 
which must be taken into account when determining the 
force necessary to separate two solids in adhesive con- 
tact. The theoretically most well defined situation is the 
contact between a ball and a flat surface. In this case, 
if at least one of the solids is elastically soft enough, the 
so called JKR theory can be applied, which predict the 
pull-off force to be [5] 

F = yTeffi?. (6) 

In this equation is the effective interfacial energy 
(usually denoted as the work of adhesion) which includes 
the influence of the surface roughness on the interfa- 
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FIG. 3: The interfacial binding energy 7 e ff (in units of the 
interfacial binding energy A7 = 71+72—712 for flat surfaces) 
as a function of the root-mean-square roughness. The red line 
is the full theory (from Fig. [2j) while the green line is the result 
for the interfacial energy assuming full contact {A — Aq) for 
all rms roughness values. 



(a) (b) 




FIG. 4: (a) An opening of crack during pull-off. (b) A local 
detachment event in front of the crack tip. Asperity contact 
region detach rapidly with elastic energy radiating inside the 
block rather than being used to break other asperity contact 
regions (schematic). 

cial binding energy. During pull-off the interfacial bond- 
breaking will not occur everywhere simultaneously, but a 
(circular) interfacial opening crack will propagate at the 
interface, see Fig. H^a). The energy to propagate an in- 
terfacial crack may be larger than the energy necessary 
just to break the interfacial bonds. Thus, the work of ad- 
hesion which enter in (6) may be larger than the effective 
interfacial energy given by (3). For example, for rubber 
materials, unless the crack moves extremely slowly, there 
may be a large energy dissipation in the rubber close to 
the crack tip caused by the viscoelasticity of the rubber, 
which may enhance the force necessary for pull-off by sev- 
eral order of magnitudes. This effect can be taken into 
account by multiplying the 7 e ff given by (3) with a factor 
l + f(v, T), where f(v, T) is caused by viscoelastic energy 
dissipation in front of the crack tip, which depends on the 
crack tip velocity v and the temperature rjjjl. l23|. And 
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even for perfectly elastic solids the work of adhesion may 
be larger than given by (3) since the separation of the 
surfaces at the crack tip may involve rapid events caused 
by local elastic instabilities close to the crack tip. Thus, 
for example, if a low asperity makes contact with the 
rubber the rubber may suddenly detach when the tensile 
stress from the approaching crack becomes high enough, 
see Fig. H^b). If the elastic energy stored in the vicinity of 
the contact area before detachment is radiated as elastic 
waves into the block, rather than used to break other as- 
perity contact regions, then this process will enhance the 
work of adhesion 7 e ff and the pull-off force. We will de- 
note processes of the type illustrated in Fig. [4fb) as local 
irreversible detachment processes. This mechanism has 
recently been suggested to enhance the work of adhesion, 
but the study in Ref. [24J is for a periodic surface profile 
(sinus corrugation) and it is not cleat to what extent the 
same effect contribute for randomly rough surface. 

We have studied the influence of rapid local detach- 
ment processes (Fig. H^b)) on the pull-off force for the 
same system as above but for a curved substrate sur- 
face. Thus we have curved the substrate surface so that 
it becomes a spherical cup with the radius of curvature 
R = 4410A, but the surface roughness is the same as for 
the nominally flat surfaces used above. We have mea- 
sured the force on the block during approach and sep- 
aration of the two solids, and used the JKR theory to 
determine the effective interfacial energy 7 e ff. This ap- 
proach includes the influence of local irreversible detach- 
ment processes (Fig. ffl^b)). As above, the two solids 
were first approached to the point where the total force 
on the block vanished, and then separated. We found 
that if the upper surface of the block moves with the 
velocity v = 5 m/s or less, the force on the block as a 
function of the (average) block-substrate separation was 
independent of v and the results presented below was 
obtained with v = 5 m/s. We also tested that the damp- 
ing which acted on the atoms (to remove elastic waves 
emitted from the interface) was so small that it had no 
influence on the pull-off force. For the curved surface 
without surface roughness we obtained from the maxi- 
mum adhesion force using the JKR formula, for both the 
approach and the separation, the same interfacial energy 
7eff = A7 ^100 meV/A 2 as determined above for the 
flat surfaces from the work done by the substrate on the 
bottom layer of the block. 

In Fig. [5]we show the interfacial binding energy (in 
units of the interfacial binding energy A7 = 71 +72 — 712 
for flat surfaces) as a function of the root-mean-square 
roughness. The red, blue and green lines are from Fig. 
[21 The pink and black curves (denoted "in" and "out") 
are the effective interfacial energy during approach and 
during separation, respectively, for the nominally curved 
(spherical cup) substrate (the curves was obtained from 
the maximum adhesion force using the JKR theory). 
Note that local irreversible detachment processes (in 
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FIG. 5: The interfacial binding energy 7 e ff (in units of the 
interfacial binding energy A7 = 71+72—712 for flat surfaces) 
as a function of the root-mean-square roughness. The red, 
blue and green lines are from Fig. [21 The pink and black 
curves (denoted "in" and "out") are the effective interfacial 
energy during approach and during separation, respectively, 
for the nominally curved (spherical cup) substrate (the curves 
was obtained from the maximum adhesion force using the 
JKR theory). The substrate surfaces are self-affine fractal 
(randomly rough) with the Hurst exponent H = 0.8. 

front of the crack tip) give a relative small increase in 
the pull-off force for randomly rough surfaces, and have 
a negligible influence on the force on the block during 
contact formation (approach). 

The rapid decrease in 7 e fj observed in the MD sim- 
ulations for small surface roughness (rms < 0.1 nm) 
may reflect atomistic binding effects. That is, the mini- 
mum energy state for the flat surfaces correspond to the 
block atoms occupying hollow sites on the substrate. For 
the rough surface other binding configurations may form, 
e.g., domain- wall structures which would tend to increase 
the elastic energy stored at the interface and hence reduce 
7eff. Such atomistic effects are, of course, not included 
in the elastic continuum model. 

Fig. [6] shows snapshot pictures of the contact region 
during pull-off for contact between the curved substrate 
with the rms roughness 3 A [(a) and (b)] and 9 A [(c) and 
(d)]. Fig. [6fa) and (c) are for the point where the total 
force on the block vanish and (b) and (d) for 0.14 ns 
later where the upper surface of the block has moved 
7 A upwards. Note that the contact region for the rms 
roughness 3 A is compact, i.e., A/Aq = 1 where A$ is the 
nominal contact area, while for the rms roughness 9 A the 
contact is incomplete. This is consistent with the theory 
prediction [see Fig. EJb)] which shows complete contact 
for the rms roughness 3 A but only partial contact for 
the rms roughness 9 A. 

Finally, we note that according to the continuum me- 
chanics theory (see, e.g., [15[), adhesive contact me- 
chanics depends on the dimensionless parameter 6 = 
E*hr ms qo/Aj and we could have plotted j^/Aj as 
a function of instead of h rms to emphasize the 
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for randomly rough surfaces local irreversible detachment 
processes (in front of the crack tip) have a relative small 
influence on the pull-off force and a negligible influence 
on the force on the block during contact formation (ap- 
proach) . 



FIG. 6: Snap shot pictures of the contact region during pull- 
off for contact between the curved substrate with the rms 
roughness 3 A [(a) and (b)] and 9 A [(c) and (d)] and block 
surface. Pictures (a) and (c) are for the point where the total 
force on the block vanish and (b) and (d) for 0.14 ns later 
where the upper surface of the block has moved 7 A upwards. 



more general or universal nature of our results. Of 
course, the result also depends on the detailed form of 
the surface roughness power spectra used. The MD 
model is atomistic and depends on an additional dimen- 
sion less number, namely the Tabor number [25[ \i = 
J R 1 / 3 A7 2 / 3 ^ _2 / 3 r " 1 , where r is an atomic bond dis- 
tance. Estimation of the Tabor number for our case 
shows that for the curved surface \i ~ 5 so we are basi- 
cally in the JKR-regime and therefore the pull-off force is 
given to a good approximation by Eq. (6) rather than the 
result of Derjaguin (valid for rigid solids) F = 27ri?7 e ff , 
or some intermediate value (e.g., DMT-theory) [26] . This 
is also confirmed by the fact that for smooth surfaces the 
A7 deduced for the fiat substrate surface is the same as 
the A7 deduced from the JKR-formula for curved sub- 
strate surface (if the Derjaguin theory would instead be 
valid, the A7 deduced using the JKR theory should be 
a factor of 4/3 times larger than deduced for the flat 
surface). 

To summarize, we have presented a molecular dynam- 
ics (MD) study of the adhesive contact between elastic 
solids with randomly rough surfaces. We have calculated 
the interfacial binding energy and compared the results 
with the predictions of a recently developed contact me- 
chanics model, which is based on continuum mechanics. 
There is good general agreement between the MD-results 
and the continuum mechanics theory, and the observed 
differences can be attributed to finite size effects and to 
the long range wall- wall interaction used in the MD simu- 
lation, in contrast to the infinite system size and contact 
interaction assumed in the theory. We have found that 
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